imgPath1 = 'D:\JC2021\IVUS menkong\dataset\标签200\';        % 标签路径
imgPath2 = 'D:\JC2021\IVUS menkong\dataset\原图200\';        % 原图路径
imgDir1  = dir([imgPath1 '*.png']);                % 遍历所有jpg格式文件
imgDir2  = dir([imgPath2 '*.png']);
pix_sum = [];
a = [];
% for I = 1:length(imgDir1)           % 遍历结构体就可以一一处理图片了
for I = 1:length(imgDir1)
    image_name = imgDir1(I).name;   % 从数组中获取文件名
    img = imread([imgPath1, imgDir1(I).name]);  % 读取标签
    Img = imread([imgPath2, imgDir2(I).name]);  % 读取原图
    img = im2bw(img);
    c = regionprops(img,'Centroid');
    centroid = cat(1,c.Centroid);  % 获取标签质心坐标
    x_label = round(centroid(1,1));  % 读取中心的x坐标，四舍五入取整，便于后面索引
    y_label = round(centroid(1,2));  % 读取中心的y坐标，四舍五入取整，便于后面索引
   
    pix = Img(y_label-190:y_label+190,x_label); % 输出上下190像素的一条，拼在一起成为一张图
    
    
    pix_sum = [pix_sum,pix];
end



imshow(pix_sum),hold on


